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low temperatures and mean number of photons, provided the damping constant is 
large enough. We investigate the quantum-classical border, the entanglement and 
K^ , decoherence of an analytically solvable model, analog to the atom-cavity system, in 

ljr\ • which the atom (field) is represented by a (driven and damped) harmonic oscillator. 

The maximum value of entanglement is shown to depend on the initial state and 
the dissipation-rate to coupling-constant ratio. While in the original model the 
atomic entropy never grows appreciably (for large dissipation rates), in our model it 
reaches a maximum before decreasing. Although both models predict small values of 
QQ ' entanglement and dissipation, for fixed times of the order of the inverse of the coupling 

^^ ■ constant and large dissipation rates, these quantities decrease faster, as a function of 

the ratio of the dissipation rate to the coupling constant, in our model. 

Ij This research was partially funded by DIB and Facultad de Ciencias, Universidad 
Nacional de Colombia. 



PACS numbers: 03.65.Ud, 03.67.Mn, 42.50.Pq, 89.70.Cf 



Damped Driven Coupled Oscillators 2 

1. Introduction 

One expects quantum theory to approach to the classical theory, for example in 
the singular limit of a vanishing Planck's constant, /i — * 0, or for large quantum 
numbers. However, dissipative systems can bring forth some surprises: for example, 
QED (quantum electrodynamics) cavity fields interacting with two-level systems, may 
exhibit classical or quantum behavior, even if the system is kept at very low temperatures 
and if the mean number of photons in the cavity is of the order of one [U, [2] , depending 
on the strength of the damping constant. Classical behavior, in this context refers 
to the unitary evolution of one of the subsystems, as if the other subsystem could be 
replaced by a classical driving. In QED cavities, the atom, which enters in one of the 
relevant Rydberg states (almost in resonance with the field sustained in the cavity), 
conserves its purity and suffers a unitary rotation inside the cavity - exactly as if it 
were controlled by a classical driving field - without entangling with the electromagnetic 
field. This unexpected behavior was analyzed in reference [1] employing several short- 
time approximations, and it was found that in the time needed to rotate the atom, its 
state remains almost pure. 

Other driven damped systems, composed by two (or more) subsystems can be 
readily identified. Indeed, in the last years there has been a fast development of quite 
different physical systems and interfaces between them, including electrodynamical 
cavities [3l H], superconducting circuits [5l |6], confined electrons [3 [H [9] and 
nanoresonators [IDl [HI [12], on which it is possible to explore genuine quantum effects 
at the level of a few excitations and/or in individual systems. For instance, the 
interaction atom-electromagnetic field is exploited in experiments with trapped ions 
[l3l[Tlj, cavity electrodynamics and ensembles of atoms interacting with coherent states 
of light [15], radiation pressure over reflective materials in experiments coupling the 
mechanical motion of nanoresonators to light [l2], and the coupling of cavities with 
different quality factors in the manufacturing of more reliable Ramsey zones [16]. In 
many of these interfaces it is possible to identify a system which couples strongly to the 
environment and another which couples weakly. For example, in the experiments of S. 
Haroche the electromagnetic field decays significantly faster [17] (or significantly slower 
[T6] ) than the atoms, the quality factor Q of the nanoresonators is much smaller than 
that of the cavity, and the newest Ramsey zones comprise two coupled cavities of quite 
different Q. Several of these systems therefore, can be modelled as coupled harmonic 
oscillators, one which can be considered dissipationless. 

In this contribution we study an exactly solvable system, composed of two 
oscillators, which permits the analysis of large times, shedding additional light on the 
classical-quantum border. Entanglement and entropy, as measured by concurrence and 
linear entropy, are used to tell "classical" from quantum effects. 
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2. The model 

The system that we consider in this manuscript comprises two oscillators of natural 
frequencies uji and uj2., coupled through an interaction which conserves the (total) 
number of excitations and whose coupling constant abruptly changes from zero to g 
at some initial time, and back to zero at some final time. We take into account that the 
second oscillator loses excitations at the rate 7, through a phenomenological Liouvillian 
of Lindblad form, corresponding to zero temperature, in the dynamical equation of 
motion [18] . Lindblad superoperators are convenient because they preserve important 
characteristics of physically realizable states, namely hermiticity, conservation of the 
trace and semi-posit ivity [19]. In order to guarantee the presence of excitations, the 
second oscillator is driven by a classical resonant field. 

The interaction can be considered to be turned on (off) in the remote past (remote 
future) if it is always present (coupled Ramsey zones or nanoresonators coupled to 
cavity fields), or can really be present for a finite time interval (for example in atoms 
travelling through cavities). The initial states of the coupled oscillators also depend 
on the experimental setup, varying from the base state of the compound system to 
a product of the steady state of the coupled damped oscillator with the state of the 
other oscillator. Since we want to make comparisons with Ramsey zones, the choices in 
the formulation of this model have been inspired by the analogy with the atom-cavity 
system, for a cavity -kept at temperatures of less than IK- whose lifetime is much 
shorter than the lifetime of Rydberg states, allowing us to ignore the Lindblad operator 
characterizing the atomic decay process. The first oscillator therefore is a cartoon of 
the atom, at least in the limit where only its first two states are significantly occupied, 
while the second oscillator corresponds to the field. 

All the ingredients detailed before can be summarily put into the Liouville-von 
Neumann equation for the density matrix p of the total system 

— = — — [i^,p] + 7(2apa^ — a^ap — pa^a) (1) 

UjL ft 

where H is the total Hamiltonian of the system and the second term of the rhs of ([1]) 
is the Lindblad superoperator which accounts for the loss of excitations of the second 
oscillator. In absence of the coupling with the first oscillator, the inverse of twice the 
dissipation rate 7 gives the mean lifetime of the second oscillator. The first two terms 
of the total Hamiltonian 

H = hiuib^b + nLU2a^a + hg {Q{t) -Q{t + T)) {a^b + dP) + ihe{e-''^°^a) - e*'^^*a), (2) 

are the free Hamiltonians of the two harmonic oscillators; the next term, which is 
modulated by the step function 0(t), is the interaction between them and the last is 
the driving. The bosonic operators of creation b (a) and annihilation b'^ {a)) of one 
excitation of the first (second) oscillator, satisfy the usual commutation relations. From 
here on we focus on the case of resonance, uoi = uj2 = ojd = ^- The interaction time 
T is left undefinite until the end of the manuscript, where we compare our results with 
those of the atom-cavity system. 
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3. Dynamical evolution 

The solution of the dynamical equation ([T]) can be written as 

pit) = V{(3{t), a{t))~p{t)V\m, oiit)), (3) 

where V{(3{t), ait)) is the two-mode displacement operator, 

V{f3{t),a{t)) = VMt))V2{a{t)) = g/^W^^-z^'W^e^Wat-a-Wa^ 

and pit) is the total density operator in the interaction picture defined by equation ([3]). 
By replacing ([3]) into ([1]), and employing the operator identities 

d ^, ^ f a*a — Q*a . „| .*^\^, n ^, sfa*a — a*a . „| 
—T>ia) = I h aa' — a a \ T>ia) = T>ia) I \- aa^ — a a 

with the dot designating the time derivative as usual, we are able to decouple the 
dynamics of the displacement operators, obtaining the following dynamical equations 
for the labels a and (3 



d I a \ I —7 — iuj 

dt \ (3 j \ —ig —iuj 




(4) 



for times between zero and T. On the other hand, the Ansatz ([3]) also provides the 
equation of motion for p(t), which turns out to be very similar to ([1]) but with the 
hamiltonian H = Hie = 0), that is, without driving. The separation provided by our 
Ansatz is also appealing from the point of view of its possible physical interpretation, 
because the effect of the driving has been singled out, and quantum (entangling and 
purity) effects are extracted from the displaced density operator pit). 

The two oscillators interact after the second oscillator reaches its stationary 
coherent state 



P2(t)=trip(t)= -e— *W-e"'-* , (5) 

as can be verified by solving (jlj) with the interaction turned off. If we want a mean 
number of excitations of the order of one then the driving amplitude must satisfy 6^7, 
and thereby the larger the dissipation is, the larger the driving is to be chosen. At zero 
time, when the oscillators begin to interact, the state of the total system is separable 
with the second oscillator state given by ([5]). The first oscillator, on the other hand, 
begins in a pure state which we choose as a linear combination of its ground and first 
excited states (again inspired on the analogy with the atom-cavity system). Thus, the 
initial state /3(0) given by 

V (0, -\ fcos(g) |0) + sin(g) |1)) (cos(g) (0| + sin(g) (1|) |0) {0\ V^ (o, -\ , (6) 

m 

corresponds to a state of the form described by equation ([3]) with /3(0) = and 
a(0) = e/7. At later times, the solution maintains the same structure, but -as can 
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be seen from the solution of (jl]) - the labels of the displacement operators evolve as 
follows 

a{t) = ee-^(^+2Ht f 1 ^^^ ^~^) ^ sin^ j ^ 

m = -^—^ + -e-^(-+2^-)* (cos {gt) + (-2/ + f ) ^^| , (8) 

where we have defined the new constant g = |y4^^ — 7^. We employ g, which also 
appears in the solution of the displaced density operator, to define three different 
regimes: underdamped {g^ > 0), critically damped {g"^ = 0) and overdamped {g'^ < 0) 
regime. It is important to notice that there is no direct connection with the quality 
factor of the damped oscillator: it is possible to have physical systems in the overdamped 
regime defined here even with relatively large quality factors, if the interaction constant 
g is much smaller than u, the frequency of the oscillators. 

The inspection of the equations ([7]) and ([8]), allows one to clearly identify the time 
scale 2/7, after which the stationary state is reached and the state of the first oscillator 
just rotates with frequency u and have a mean number of excitations equal to e^/g"^- 
The doubling of the damping time of the second oscillator, from I/7 in the absence 
of interaction to 2/7, in the underdamped regime, can be seen as an instance of the 
shelving effect [20]. The first oscillator, which in absence of interaction, suffers no 
damping, it is now driven and damped. It can be thought that the excitations remain 
half of the time on each oscillator, and that they decay with a damping constant 7, 
thereby leading to an effective damping constant of 7/2. An interesting feature of the 
solution is that the displacement of the second oscillator goes to zero, in the stationary 
state. In the stationary state, the first oscillator evolves as if it were driven by a classical 
field —ihe exp{—iujt) and damped with a damping rate g, without any interaction with 
a second oscillator. More generally speaking, we remark that from the point of view of 
the first oscillator, the evolution of its displacement operator happens as if there were 
damping but no coupling, and the driving were of the form hg{(3 — ia), or, in terms of 
the parameters of the problem, 

F{t) = -thee-''^' - znee-(^/'+*")* ((~ ~ "^^^^^^^ + ^^ +91-1 sin(^^)'\ , (9) 

This behavior is particularly relevant in the following extreme case, whose complete 
solution depends only on the displacement operators. If the initial state of the first 
oscillator is the ground state then p does not evolve in time, i. e. it remains in the state 
|00), and the total pure and separable joint state is 

Pit) = mt)) {m\ ® Ht)) Ht)\ . (10) 

Even in the more general case considered here, corresponding to the initial state 
([6]), the solution of p{t) possesses only a few non-zero elements. If we write the total 
density operator as 

m= E Pn^l^i^2)(jU2|, (11) 



*1«2J1J2 
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we can arrange the elements corresponding to zero and one excitations in each oscillator, 
as the two-qubit density matrix 

f pW pW pW o\ 



V 



(12) 



P'o'iit) P'olit) pllit) 
pToit) Pirn pio(t) 

y 

If we measure time in units of g by defining t' = gt we have only two free parameters 
r = - and Q = -. The nonvanishing elements of the density matrix, written in the 
underdamped case (|r| < 2), are given by (hermiticity of the density operator yields the 
missing non-zero elements) 



P'o'oit') 



1 — sin 



-rt' 



, (A - r2 cos(V4 - r2t') r sin(V4 - r^t') 



4-P 



V4^^T2 



Poi(tj - ism (6^)6 ^_-p2 

-10/ >^ _ ■ 2M^ ~Tt' I (2 - r') cos {VI^^T^t') +2 r sin {^/I^^T^t') 



pl^it') = sin^(0)e- 



p^^(f) =isin(2^)e 



iUf 



I rt' 



4_r2 

sin (y/I^rpf) 
v/4^^T2 



V4^^ 



~oo/^/^ _ sin(2^^if^i._rv 



p1[;(t' 



cos 



yl^^p- 



t'\ rsin(v/4^^P|)' 



pi'i(t') = 2^sin'(^)e 



^ri 



2/ V4^^T2 

,sin(v/4^^P|) /rsin(v/4^^n|) 



V4^^ 



VI^^ 



cos 



VI^^P- 



The expressions of the elements of the density matrix in the critically damped case 
r = 2 and in the overdamped case F > 2 can be obtained from those given in the text 
for the underdamped case F < 2. 



4. Entanglement 

Although quantities like quantum discord [21] have been proposed to extract the 
quantum content of correlations between two systems, we presently quantify the 
quantum correlations between both oscillators employing a measure of entanglement. 
Due to the dynamics of the system, and the initial states chosen, the whole system 
behaves as a couple of qubits and therefore its entanglement can be measured by 
Wootters' concurrence p^. One of the most important characteristics of the form 
of the solution given by ([3]) is that concurrence, as well as linear entropy, depend only 
on the displaced density operator pit'). In our case the concurrence reduces to 



C{t') 



p%{trpf,it') + \ p'^it'yp'^it' 



Vp?o(i')pl?(t') = 2|p^?(t')|, 



p%{trp}^it') - jpgi(t')pig(t' 



(13) 
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where the positivity and heriniticity of the density matrix were used. The exphcit 
expressions for the concurrence in the underdamped (UD), critically damped (CD) and 
over damped (OD) regimes are 



CuD{t')=^sin\e)e 



-rt' 



sm 



Vi^Ta t' 



V4-r2 

CcD{t') = 2sin^{e)e-^W\t' -l\, 

sinh 



Fsin 



v^i^Ta t' 



VI^^ 



COS 



v/I^^Pt' 



(14) 



CoD(t')=4sin2(^)e-r*'- 



v^rani t' 
2 



vT234 



Fsinh 



vr^^ t' 



vT2 



cosh 



Vn"^t' 



All the dependence on the initial state is contained on the squared norm of the coefficient 
of the state |1) of the displaced density operator. In all regimes the concurrence vanishes 
at zero time, because the initial state considered is separable. However, while in the 
underdamped case the concurrence vanishes periodically (see equation (IT^ below), in 
the other two cases it crosses zero once (t > 0) and reaches zero assymptotically as time 
grows. This shows a markedly different qualitative behavior (see figures [1] and [2]) . 
In the underdamped regime the zeroes of the concurrence are found at times 

2n7r 27rn + 2arccos (^) 

Tin = , and T2n = , -, (15) 

where n is a non-negative integer. In this contribution, the inverse sine and cosine 
functions are chosen to take values in the interval [0, 7r/2]. The time tiq corresponds to 
the initial state. The sequence of concurrence zeroes is thereby = rio < r2o < th < 
T21 . . . As the critical damping is approached, the time tu is pushed towards infinity, 
while T20 approaches the finite time 2/r (see figure [2]). For the initial states considered 
in this manuscript we do not observe the sudden death of the entanglement since the 
concurrence is zero only for isolated instants of time. 

If one writes the concurrence in the underdamped regime in the alternative form 

sin^^ 



Cunit') 



-rt' 



2(1-FV4) 



sm 






arcsm — 




it is easy to verify that at the times T±n, given by 



r± 



{2n + l)7r ± arccos 



F^ 



the concurrence reaches the local maxima 

T{{2n 



2 arcsin 



>0, 



n 



0,1, 



(16) 



(17) 



a 



±n 



sin'^ 




exp 



l)7r ± arccos ( 



— 2 arcsin 



1-? 




We observe these maxima to lie on the curves sin^ 6'i^-i- exp(— Ft), where the constants 

K± = Jl + ^± J satisfy the inequalities V2 - 1 < K. < 1 < K+ < V2 + 1. Maxima 
of concurrence depend on both the initial state and the value of the rescaled damping 
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constant, and reach the maximum available value of one only in the non-dissipative case 
for a particular initial state. In order to have negligible values of concurrence (except for 
small time intervals around the zeroes of concurrence) it is necessary to have times much 
larger than I/7. From the point of view of classical-like behavior, the most favourable 
scenario corresponds to zero or almost zero concurrence, which are obtained for short 
time intervals around ri„, r2„ and for large values of time. 

In the overdamped regime, the concurrence presents two maxima, r_ and r+ > t^ 

2arccosh(r/2) ± arccosh(rV4) ,^ „, 

'T± = , = , (18) 

both of which go to zero as the rescaled dissipation rate grows, r+ —>■ 41n(r)/r and 
r_ —>■ ln(2)/r (see figure [2]). The function arccosh(a;) is chosen as to return nonnegative 
values for x > 1. Since the global maximum of concurrence, which corresponds to the 
later time, scales like l/(2r) for large values of F, in the highly overdamped regime 
quantum correlations are not developed at any time. 




0.1? 
0.10 
O.O.s 



A 










\ 




(b) 






\ 










\ 


^ 


\ 






\ 


1 


\ 







2 3 4 5 

t 




Figure 1. Concurrence as a function of time and the rescaled damping constant in 
the (a) underdamped, (b) critically damped, and (c) overdamped case. 

The behavior of concurrence in the different regimes is shown in figure [H It is 
apparent that small values of concurrence are obtained for very small times and for large 
times in the underdamped case and for all times for a highly overdamped oscillator. In 
figure [2] we depict the times at which concurrence attains a maximum, and the maximum 
values of concurrence, as a function of the rescaled damping constant. One can see how 
the first two times of maximum concurrence go to zero, while the other times diverge, 
as the critically damped regime is reached. The first two maxima of concurrence vanish 
more slowly than the rest of maxima, which hit zero at F = 2. 

5. Entropy 



The entropy is analized employing the linear entropy of the first oscillator, the system of 
interest. As remarked before, the first oscillator behaves like a two- level system, where 
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(b) 



Figure 2. (a) Concurrence local niELxinia at times !+(_)„ in black (gray) color for 
n = full line, n = 1 dashed line and n = 2 dashed-dot line and (b) times of maximum 
values of concurrence, as a function of the rescaled damping constant. 



the maximum value of the hnear entropy, 0.5, is obtained when the population of each of 
the two states is one half. The type of "classical" behavior which allows the interaction 
with Ramsey zones to be modelled like a classical driving force occurs when the linear 
entropy is very small, and hence the state of the first oscillator is (almost) pure and 
uncorrelated with the state of the second oscillator. The linear entropy for the first 
oscillator can be computed as 

6i{t') = l- tri (trs p trs p) = 1 - tri (trs p tra p) = 2 det (trs p) , (19) 

where the last equality holds for two-level systems. In equation (IT9l) the density operator 
of the first oscillator is assumed to be represented by a 2x2 matrix. Employing the 
expressions we have found for the elements of p we obtain 

6i{t') = 2 sin^ e x{t') (1 - x{t')) , (20) 

where x, in the underdamped regime, is given by 

(^l-V^/A t - arccos(r/2)' 



— r^ • 2 
e sm 



xuoit) 



(21) 



1 - r74 

Surprisingly, as in the case of concurrence, the influence of the initial state factors out 
in the expression of the linear entropy of the first oscillator, which turns out to be 
proportional to the square of the population of |1) in the initial displaced operator. 
As it is well known, in the limit of zero dissipation, the linear entropy of the reduced 
density matrix is equal to one fourth of the square of the concurrence. At times r2„ 
(see eq. (ll5l) ). when both concurrence and linear entropy vanish, the total state of the 
system is separable, p{gr2n) = \P{gT2n)) (/3(5'T2n)| ® P2(fl'T2n); that is, from the point 
of view of the first oscillator the evolution is unitary like. Since the linear entropy 
begins at zero, because the initial state is pure and separable, there is a maximum in 
the interval (0,r2o), which turns out to give a linear entropy of exactly 0.5 (we treat 
the case sin 6' = 1, because — due to the scaling property discussed before — a simple 
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multiplication by sin^^ gives the result for other cases). Indeed, as the function x{t) 
changes continuously from x{t = 0) = 1 to x{t = T20) = 0, it crosses 0.5 at some time 
Tso in between, giving the maximum value possible of the linear entropy. Although the 
exact value of t^q can be obtained only numerically, good analytical approximations can 
be readily obtained. For example, t^q ~ tt/{4 + Ag + 2g'^), gives an error smaller than 
0.5%. 

For small values of the rescaled damping constant, F ^ 0.237, there are several 
solutions to the equation x(t) = 0.5 in the interval (0,log(2)/F), which give absolute 
maxima of the linear entropy, while the times 

r,„^^-^E^£M±^, n = 0,l,2,... (22) 

correspond to local minima. In the interval (log(2)/F, cxd) the times r^n give local 
maxima. All of the local maxima and minima given by eq. fl22|) belong to the curve 
2e~'"*(l — e~^*). The large time behavior of the local maxima of linear entropy and 
concurrence is, thereby, of the same form constantx exp(— Ft'). For values of F > 0.237 
all times r4„ give local maxima. The maxima of concurrence and linear entropy coincide 
only in the weakly damped case, because concurrence and linear entropy are not 
independent for pure bipartite states. 

At times Tin (see eq. flT^ ). where the total state p{grin) = Pi{gTin) ® 
\a{gTin)) {a{gTin)\, is separable, the reduced state of the first oscillator is mixed. The 
linear entropy is small for short (ri„ ^ log(2)/F) and large (ri„ ^ log(2)/F) times. 

In the overdamped regime the function x{t), which appers on the expression for 
linear entropy ( l20l) and is given by 

e-r*' sinh2(v/FV4-l t' - arccosh(F/2)) 

^OD[t ) - fVl^l ' ^^ 

begins at one for t' = 0, and goes down to zero for large values of time. The time at 
which it crosses one half can be calculated to be ro.5 ~ 0.16557 < 1/6 for F = 2 and 
for large values of F it goes as tq.s ~ In2/(2F). It is easy to find interpolating functions 
with small error for the time of crossing, 

/ \ 

1 



pI'o 



t' 



\ ^ "*" In2 ^^'^^ ( arccosh(^) tanh I ^-g — — 1 1 1 



0.5(1 + A) 



where |A| < 2.5%. It is interesting to notice that for large values of the damping this 
time (ln(2)/(2F)) is half the time needed to obtain the maximum value of concurrence, 
and that, at the later time, the linear entropy is 3/4 of the maximum value of entropy, a 
relatively large value. The state of the first oscillator always becomes maximally mixed 
before becoming pure again, no matter how large the value of the damping. We show 
the behavior of linear entropy in figure [31 In the underdamped regime there are infinite 
maxima and minima, while for critical damping and for the overdamped regime there 
are only two maxima. The first maximum always corresponds to a linear entropy of one 
half. 
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b 



12 3 4 




Figure 3. Linear entropy of the first oscillator as a function of time and the rescaled 
damping constant in the (a) underdamped, (b) critically damped, and (c) overdamped 



case. 



6. Conclusions 

In the present contribution we have shown that the classical quantum border in this 
model depends mainly on the initial state and on damping constant to interaction 
coupling ratio, and that quantum effects, characteristic of the underdamped regime, 
can be seen in the other regimes for small times. In order to make connection with 
Ramsey zones we remember in that physical system uj k, 10^'^ Hz, Q k, 10^, g ~ 10^ 
Hz and Tr k. 10~^ s, which was chosen as to produce 7r/2 pulse, that is a pulse that 
can rotate the state of the two-level system, as represented in a Bloch sphere, by an 
angle 7r/2. These numbers place the system into the highly overdamped (regime because 
r = uj/{Qg) ~ 10^ S> 2) and give a rescaled evolution time of the order of gT pa 10~^. 
Here we use the same values of uj, 7 and g, and an evolution time of order 1/g. Indeed, 
the hamiltonian hub^b + hg (e(t) - e(t + T)) (aoe"*'^*&"^ + aje^^'fo), with ||ao|| ~ 1 — 
which would model the interaction of the first oscillator with a classical driving field of 
an average number of excitations of the order of one — has a characteristic time 1/g, 
corresponding to T' ^ 1. 

The dynamical behavior of the linear entropy obtained here, is quite different 
from that of ref. [1]: there the linear entropy was never large for the relevant time 
interval, here it grows to the maximum possible for a two-level system, and then goes 
to zero very quickly. Therefore, in this model dissipation produces relaxation also, and 
a description obviating the second oscillator still needs a dissipation process. Although, 
at the evolution time T, both models predict a small atomic entropy, in Ramsey zones 
it decreases as 5i(Tr' ~ 0.1) ~ 4/r, while in the present model it goes to zero as 
Si(T' ^ 1) oc 1/r^. Qualitative and quantitative differences notwithstanding, at the 
evolution time the linear entropy is very small, in both cases, due to the smallness of 
the ratio g/'j. As remarked before the quality factor of the damped oscillator does 
not appear directly in either case; it can be perfectly possible to have a very weakly 
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damped oscillator and a highly overdamped interaction. However, as the first oscillator 
quality factor is improved, the damping constant will eventually be comparable with 
the interaction constant, and there will be considerable entanglement between both 
oscillators. For the same physical system if the damping rate can be changed then 
classical or quantum behavior can be obtained. 
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